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We investigate the properties of a cosmological scenario which undergoes a gravitational phase 
transition at late times. In this scenario, the Universe evolves according to general relativity in 
the standard, hot Big Bang picture until a redshift z < 1. Non-perturbative phenomena associated 
with a minimally-coupled scalar field catalyzes a transition, whereby an order parameter consisting 
of curvature quantities such as R 2 , R a tR ab , R a bcdR abcd acquires a constant expectation value. The 
ensuing cosmic acceleration appears driven by a dark-energy component with an equation-of-state 
w < — 1. We evaluate the constraints from type la supernovae, the cosmic microwave background, 
and other cosmological observations. We find that a range of models making a sharp transition to 
cosmic acceleration are consistent with observations. 
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Introduction. The observational evidence for a low den- 
sity, spatially- flat, accelerating Universe poses a severe 
challenge to theoretical physics. Solutions which have 
been proposed include a cosmological constant, an ultra- 
light scalar field, and modifications to Einstein's gravity. 
An approach which unifies these different viewpoints can 
be found in Sakharov's description of the gravitational 
physics of the vacuum 0. According to his proposal, 
quantum effects of the particles and fields present in the 
Universe give rise to the cosmological constant and grav- 
itation itself. Although this approach has not solved the 
problem of the cosmological constant, it has pointed the 
way towards a novel model for the dark energy. 

Quantum effects of an ultra-light, minimally-coupled 
scalar field have been proposed to account for the dark 
energy phenomena in the vacuum metamorphosis sce- 
nario [2 El [3 IE Effectively, the Ricci scalar cur- 
vature serves as an order parameter, marking a gravita- 
tional transition when it drops to the value xm 1 . Here, 
m ~ 10 -33 eV is the mass of the scalar field and % is 
a numerical constant of order unity which is set by the 
theory. For most of the history of the Universe, up until 
z ~ 1, R > xm 2 and the vacuum stress-energy is neglig- 
ble. Since the local value of the Ricci scalar exceeds this 
critical value in the vicinity of galaxies today, we see no 
vacuum energy nearby. On the largest scales, however, 
the vacuum term comes to dominate the cosmic evolu- 
tion as R — > vm 2 - A detailed analysis shows that feed- 
back from the vacuum prevents R from evolving past the 
critical value, in a sort of gravitational Lenz' law which 
maintains R — x m2 Ifll- The ensuing, large-scale cosmic 
expansion begins to accelerate, driven by a dark energy 
component with an equation-of-state w < — 1 which ap- 
proaches — 1 in the future. 

Despite the sharp change in the character of gravita- 
tion, we are still justified to use basic cosmological tools 



such as the luminosity distance. Although the transition 
changes the relationship between curvature and matter, 
we still have a metric theory and a complete description 
of the evolution of the expansion scale factor. 

The vacuum metamorphosis scenario is distinct from 
scalar-tensor theories of gravity or high-energy 
physics inspired modifications of the gravitational action 
whereby the Einstein-Hilbert action is replaced by a func- 
tion of the Ricci scalar, — l6irG£ g = R — > f(R). (See 
Ref. [8j for a recent example.) In our case, the gravita- 
tional action is modified by vacuum polarization effects 
of the matter content in the theory. However, the modifi- 
cations include contributions from R, R a b, Rabcd so there 
is not a simple way to re-express the model in terms of 
an equivalent, non-minimally-coupled scalar field. (See 
Refs. 13, [l0| for a discussion of the equivalence between 
higher-order gravity theories and scalar-tensor gravity.) 
Furthermore, the scenario we consider is different from 
other dark energy models based on curved-space quan- 
tum field theory effects 0, E3 > in that the cosmic ac- 
celeration arises from non-linear, non-perturbative phe- 
nomena. In fact, the theory we examine in this article 
is similar in spirit to the gravitational transition inves- 
tigated by Tkachev yj|]. The distinction here is in the 
form of the vacuum stress-energy, which is based on the 
novel results by one of us and collaborators 0, IH 0, 0, 0] • 
The one-loop effective action, which is complete for a free 
massive scalar field in curved spacetime, can be viewed 
perturbatively in terms of Feynman diagrams. Classical 
gravitons attach to the vacuum loops of the scalar field. 
With increasing number of external gravitons, these di- 
agrams require counterterms that give rise to renormal- 
ization of the cosmological constant, Newton's constant, 
and several other constants that appear in terms that are 
of second order in the Riemann tensor. In the limit of a 
massless field, the conformal trace anomaly arises from fi- 
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nite parts of the second order terms and does not depend 
on the values of the renormalized constants. The terms in 
the perturbative expansion that are of still higher order 
in the Riemann tensor have no infinities and appear to be 
relevant only at Planckian scales. However, summation 
of an infinite subset of those terms shows that they may 
give rise to a large nonperturbative effect when the Ricci 
scalar curvature approaches a specific value proportional 
to the square of the scalar particle's mass. Alternatively 
the functional integral over fluctuations of the scalar field 
in the one-loop effective action can be performed and the 
result can be expressed in terms of the exact heat kernel 
of the scalar field equation. By studying known exact 
solutions for the heat kernel, one can infer a plausible 
asymptotic form of the heat kernel and show that it leads 
to the same type of large nonperturbative effect when the 
Ricci scalar curvature approaches a value, xto 2 , propor- 
tional to the square of the particle's mass. This latter 
approach does not rely on a subset of Feynman diagrams 
and is thus fully nonperturbative. Like the conformal 
trace anomaly, this effect is independent of the values of 
renormalized constants. If the calculation is valid, vac- 
uum metamorphosis must occur if the universe contains 
a light scalar field with x positive. Since a complete de- 
scription of the stress-energy tensor source for Einstein's 
equations must include all matter fields as well as vac- 
uum stress-energy, the cosmological phenomena resulting 
from the vacuum stress-energy is inevitable given the ex- 
istence of such a scalar field. 

The original vacuum metamorphosis model is tightly 
constrained, although not eliminated, by observations 
[Til H5I ] . We find the basic behavior to be of sufficient 
interest to justify further investigation. In this article we 
speculate that other curvature quantities such as R a bR ah 
or R a b c dR abcd can serve a similar role as order parameters 
for a gravitational transition. In the following we exam- 
ine the cosmic evolution resulting from such gravitational 
transitions. We evaluate the observational constraints 
based on type la supernovae, the cosmic microwave back- 
ground, the Hubble constant, and large scale structure. 
We stop short of making a full perturbation analysis — 
the detailed equations require a lengthy investigation — 
and we regard this as a first cut at a family of models 
extending the vacuum metamorphosis scenario. 

Gravitational Transition. In the vacuum metamorphosis 
scenario, non-perturbative effects of a light scalar field 
lead to a gravitational transition when the Ricci scalar 



reaches the level R 



Here we absorb the dimension- 



less parameter x into the mass fh. Before the transition, 
the cosmic evolution is determined by the standard FRW 
equation 



3H 2 = %ttG P „ 



a < a* 



(1) 



where a* marks the time of the transition, and p m rep- 
resents all matter and radiation. There is no need to 



include the vacuum energy density, p v , since it is negli- 
gible at these early times. After the transition, however, 
the cosmic evolution is given by 



R = 6(H + 2H 
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a > a* 



(2) 



Notably, the subsequent evolution of matter and radia- 
tion after the transition have no influence on the expan- 
sion rate. The vacuum energy does not merely contribute 
to the cosmic energy density driving the expansion, as for 
most dark energy models. Rather, it changes the form 
of gravity on cosmological scales and completely deter- 
mines the expansion. After solving J5J) for H(a) or a(t), 
however, we can still use the standard FRW equation 



3H 2 = 8ttG (p m +Pv) 



(3) 



together with the evolution of the matter density param- 
eter 



(4) 



to deduce the properties of an equivalent dark energy 
that has the equation-of-state, w = p v /p v : 



1+^H/H 2 ) /(I -<>,„(„!). (.-,) 



The inferred equation-of-state is super-negative, with 
w < — 1 which approaches w — ► — 1 in the future. For 
this model, there is a single free parameter, fh. So, after 
choosing a present day Hubble constant H , the value of 
fh determines the matter density: 
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The Hubble constant at the time of the transition is 
fh/y/3, so 



H(a) 
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gives the analytic solution to the evolution equation J2J. 
Note that the total energy density and pressure is always 
positive: p + p = m 2 (a*/a) 4 /87rG. There is no "big rip" 
in this scenario 16j , as the late-time evolution resembles 
a comological constant-dominated universe. 

We propose to extend this gravitational transition 
model, using R a bR ab or R a bcdR abcd as the order parame- 
ter in place of R. In either case, the evolution before the 
transition is determined by the standard FRW equation 
||TJ . After the transition, from the constancy of the order 
parameter we have 



H 2 + N(HH 2 + H 4 ) 



a > a* 



(8) 
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where a — 12 and N — 2, 3 for the Riemann and Ricci 
tensors, and a = 36, N — 4 for the Ricci scalar. Linear 
combinations of these curvature terms which lead to cos- 
mic acceleration, a > 0, correspond to 3/2 < N < 9/2. 
We can solve these equations numerically to determine 
the subsequent evolution. 

The cosmic evolution under various gravitational tran- 
sition scenarios is shown in Figures 1112131 The hallmark 
of these models is a rapidly-evolving equation-of-state 
with w < — 1. We can see that by decreasing AT, or 
increasing the rank of the curvature quantities used to 
form the order parameter, the strength of the vacuum 
energy increases. The total p + p drops faster in time, 
and the onset of cosmic acceleration is sharper. 




-3.0 -2.5 -2.0 -1.5 -1.0 



equation-of-state: w 

FIG. 1: The evolution of the gravitational transition model 
in the w — w/H phase plane is shown. For a given N, the 
model trajectory follows the curve indicated with Q m decreas- 
ing as the phase variables evolve from the starting point at 
(w, w/H) = (—00,00) and head towards (—1,0). The values 
corresponding to Sl m = 0.3, 0.4 are indicated by the small 
squares and circles, respectively. In the case f2 m = 0.4 today, 
(w,w/H) = (-1.55,2.4), (-1.45,1.6), (-1.35,1.0), with the 
transition occuring at a,/ao = 0.7, 0.6, 0.5 for N = 2, 3, 4 
respectively. 

Another possibility that comes to mind as an order 
parameter is the Gauss-Bonnet invariant, 

Q = R 2 — 4:R a bR ab + R a bcdR abcd ■ (9) 

However, the deceleration parameter q = — aa/et 2 is pro- 
portional to the Gauss-Bonnet invariant, q = —Q/2AH 4 . 
This means the curvature cannot evolve down to some 
constant Q to signal a gravitational transition to cosmic 
acceleration: Q < in a matter-dominated universe, but 
Q > is required for acceleration. 

Constraints. We now consider the observational con- 
straints on the cosmic evolution resulting from the fam- 
ily of generalized vacuum metamorphosis models. These 
constraints are due to: 1) distance - redshift relationship 
using type la supernovae (SNe) 0,0]; 2) the distance 
to the cosmic microwave background (CMB) last scatter- 
ing surface and the density tt m h 2 implied by the WMAP 




FIG. 2: The evolution of the equation-of-state w as a function 
of the scale factor a is shown. All models have £l m — 0.4 
today. For decreasing N, the effective gravitational repulsion 
of the dark energy as measured by w increases. The N = 4 
model is consistent with all observations, the iV = 3 is on the 
border, and N — 2 is excluded. 
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FIG. 3: The evolution of the density parameter Q m as a func- 
tion of the scale factor a is shown. All models have Q m = 0.4 
today. The matter density drops suddenly at the onset of the 
transition. The sharpness of the drop grows with decreasing 
N. The N = 4 model is consistent with all observations, the 
N = 3 is on the border, and N = 2 is excluded. 

measurements Il9j: 3) the Hubble constant based on the 
HST key project [2fJ; 4) the mass power spectrum shape 
parameter T = Q m h |2lll22| . 

• For the SNe we use the 156 supernova, "gold" data 
set of Riess et al and the 54 supernova set of 
Knop et al We make a simple y 2 test to de- 
termine the 2cr region allowed for each data set. 

• For the CMB we exploit the geometric degener- 
acy of CDM-family models with identical primor- 
dial perturbation spectra, matter content at last 
scattering, and comoving distance to the surface 
of last scattering [2^, I24I Eif . This means that 
there is a family of dark energy models with dif- 
ferent equation-of-state histories and dark energy 
abundances but essentially identical CMB spec- 
tra [2^|. Hence, our best-fitting models are those 
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which are degenerate with the best-fitting ACDM 
models obtained by WMAP [l^. We have taken 
the 2a ACDM models based on a 5-parameter fit 
£l c h 2 i Qbh 2 , h, n* obtained from a combina- 
tion of earlier work [27|. new calculations using 
CMBfast H3, and CMBfit [H). We also use the 
Big Bang Nucelosynthesis prior Vi^ti 1 = 0.02 ± 
0.002(95%) [13 • Ultimately, these points in the 
ACDM Q rn — Hq plane are mapped to points in 
the VCDM Q m — Ho plane by fixing the quantities 
fl m h 2 and the luminosity distance to the last scat- 
tering surface. This procedure overlooks differences 
in the large-angle anisotropy pattern, which would 
require a full treatment of the cosmological per- 
turbations after the transition to model accurately. 
However, the x 2 fit between the theoretical model 
and experiment is dominated by the small-angle 
anisotropy pattern so these differences should be 
small. (A preliminary treatment of the large-angle 
CMB anisotropy suggests the VCDM constraint re- 
gion shifts to higher VL m by < 10% |l5j.) 

• For the Hubble constant, we require that the 
value of H falls within 2a of HST's result 72 ± 
8 km/s/Mpc (la) H^. 

• For the shape parameter, the chief obstacle in ap- 
plying the 2dF and SDSS results is the difference 
in the rate of growth of linear density perturba- 
tions between our models and a ACDM cosmol- 
ogy. We expect linear perturbation growth to 
cease at the transition, although we postpone a 
detailed treatment for a future investigation. The 
current la bound on T has a ~ 20% uncertainty 
which should be comparable to the differences re- 
sulting from the perturbation growth. For safe 
margin, we require that our best-fit model has a 
shape parameter which falls within 3a of the bound 
T = 0.2 ±0.03 (la) jUi^. 

The results are shown in Figure 21 We see that the 
gravitational transition models require a lower matter 
density than ACDM to satisfy the CMB constraints, and 
a higher matter density to satisfy the SN constraints. 
The \ 2 fit to the CMB data for the best-fit model is iden- 
tical to that for the best-fit ACDM model, due to the ge- 
ometric degeneracy. The fit to the SN data is marginally 
better than ACDM { x 2 = 175 for 156 SNe 0; X 2 = 59 
for 54 SNe 18]). For N = 4, the original VCDM model, 
there are viable models near Vl m = 0.4, Hq = 70 whereby 
the transition redshift is z* ~ 1 0. However, as N de- 
creases, the tension between CMB and SN data builds. 
Hence, sharpening the gravitational transition relative to 
VCDM conflicts with observation. As N approaches 9/2 
it turns out that the gap between the SN and CMB pa- 
rameter ranges decreases. But in this limit, for fixed mat- 
ter density today, the transition occurs at earlier and ear- 
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FIG. 4: The observational constraints on the Q m — Ho param- 
eter space for the N — 2, 3, 4 gravitational transition models 
are shown. The red-filled contour is the WMAP CMB 2a re- 
gion, obtained from the ACDM 2a region, given by the black 
dashed contour, by using the geometric degeneracy. The line 
connecting the black dots shows the mapping of the best-fit 
ACDM model to the gravitational transition model. The SN 
2a regions are shown by the wide, pale blue (Knop et al 2004) 
and narrow, pale grey (Riess et al 2004) bands. The shape pa- 
rameter F 3<r region is shown by the yellow swath. The Hub- 
ble constant 2a region is the horizontal green band. For ref- 
erence, gravitational transition models with age 13.5 Gyrs lie 
along the thin blue line. There is significant overlap amongst 
all model constraints for the cases N = 4 for the R 2 tran- 
sition. Sharper transitions, corresponding to N < 4, are in 
conflict with some or all of the observations. 



lier times, at ^> 1. Since linear perturbation growth is 
expected to slow or cease at the transition, such an early 
transition can be expected to wreck the basic structure 
formation scenario. 

There are some straightforward ways to ease the con- 
flict between these gravitational transition models and 
observations: 1) allow for a small amount of spatial cur- 
vature; 2) include the energy of the scalar field itself. 

Adding negative spatial curvature improves the fit to 
the SN data while lowering the matter density. However, 
the viable CMB region shifts to an even lower range of 
Q m , and so the two major constraints remain in con- 
flict. Adding positive spatial curvature, thereby closing 
the Universe, brings the 2a SN and CMB regions into 
closer agreement. As illustrated in Fig.O the addition of 
a small amount of positive spatial curvature, flk — —0.05 
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FIG. 5: Same as Fig. fj] but with positive spatial curvature. 
In all cases = —0.05 and fl m = 1.05 at the time of the 
transition. The agreement with CMB and SN constraints 
improves significantly, particularly for the N=4 case as the 
CMB best-fit point falls inside the 2a contour for all other 
constraints. 



at the transition, leads to concordance. The quality of 
the improvement is best for N = 4, whereby Q m ~ 0.4 
and Hq 60 km/s/Mpc lies within the 2a contour for all 
the constraints. Yet, we view the addition of spatial cur- 
vature as extraneous. It is unrelated to the mechanism of 
vacuum metamorphosis or gravitational transition, and 
there is no direct observational evidence which requires 
it. Hence, we do not pursue spatial curvature any further. 

Including the energy of the scalar field itself can also 
bring the SN and CMB constraint regions into better 
agreement. Here it is necessary to explain that the vac- 
uum stress-energy which gives rise to the gravitational 
transition consists of the gravitational corrections to the 
scalar field stress-energy in the state in which the vac- 
uum expectation value of the scalar field is zero. A 
non-zero vacuum expectation value would contribute the 
same stress-energy as a free classical scalar field of mass 
m. 

This classical scalar field introduces new degrees of 
freedom into the model. The mass of the scalar field, 
with potential V = \m 2 (j) 2 , is related by m 2 = m 2 /\ to 
the parameter m that determines the cosmological evo- 
lution after the transition. Using the value x = 6 sug- 
gested by earlier versions of the VCDM model 0, then 



strongly damped by the cosmic expansion up to the time 
of the transition, so that to first approximation it will 
look like a cosmological constant. In order to contribute 
a fraction 1 — f2 m (a*) of the total cosmic energy den- 
sity at the transition time, the field amplitude must be 
(f> = ^/3(1 — £! m (a*))/27rA/p where Mp is the Planck 
mass. This is a potential sticking point for the model, 
much as for scalar field quintessence, since such a high 
field amplitude should be susceptible to quantum grav- 
itational effects which, for example, induce couplings to 
all other matter fields. However, the absence of non- 
gravitational interactions for 0, together with its very 
small mass m, appear to suppress such couplings. Also, 
one may question whether it is necessary or economical 
to introduce a second form of dark energy. This reduces 
to the question of the vacuum expectation value of the 
field, which must be determined by the initial conditions 
or perhaps post-inflationary physics. Whether or not the 
vacuum expectation value is large, for the quantum scalar 
field that we have been considering, the quantum effects 
that lead to vacuum metamorphosis are inevitable. 

Including an effective cosmological constant before the 
transition, the evolution equations for the Hubble expan- 
sion rate or the scale factor are the same as (JBJ, but with 
different initial conditions at the transition. For the R 2 
case with N = 4, the analytic solution is 



H(a) 



2^3 



3fl m ((2*) 
2 



i 



1/2 



IT 



, (io) 

(11) 



The equation-of-state is still obtained from J^J. But 
whereas w is undefined at the transition in the absence of 
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scale factor: a/a 
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m /6. The evolution of such a light field will be 



FIG. 6: The evolution of the equation-of-state w as a function 
of the scale factor a is shown for the N = 4 model including 
the classical scalar field. In all cases, w starts at —1 at the 
transition, drops below, and then ultimately evolves back to 
— 1 in the asymptotic future. For Q m = 0.3 today, only the 
top curve, with Q m — 0.67 at the transition, produces a viable 
cosmology satisfying all observational constraints. 
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FIG. 7: Same as Fig. 0] but the classical mode of the scalar 
field contributes an effective cosmological constant before 
the gravitational transition. At the transition, Q m (a t ) — 
2/3, Qa(o,*) ~ 1/3. The result is that the transition to ac- 
celeration is softened, and the tension between the CMB and 
SN constraints is relaxed. Models with N — 2 — 4 are viable. 



the classical scalar field, now u>(a*) = —1 as illustrated 
in Figure Also, note that f2 m (a*) > f2 m (ao). In the 
N = 4 model with f2 m (a*) = 2/3 and fl m {ao) — 0.3, the 
transition occurs at a»/ao = 0.64 and the equation-of- 
state is w = —1.15 today. As Q m (a*) is lowered, the onset 
of cosmic acceleration is earlier but the transition is later, 
for fixed Q m (ao). In the limit that O m (a*) = Q m (ao) the 
transition occurs at z* = and the model is equivalent 
to ACDM. 

We have evaluated observational constraints for the 
R 2 , R a bR ab and R a bcdR abcd models with a cosmological 
constant. As more A is added before the transition, the 
gap between the CMB and SN constraints reduces. In 
Figure [7| we show the constraints resulting for the case 
n m (a«) = 2/3. The quality of the x 2 St to the CMB 
and SN data are approximately the same as for the cases 
without the addition of the scalar field contribution. In 
all cases shown there are viable models near VL m = 0.3 
and Hq — 70, providing good motivation for further in- 
vestigation of these models. 

Cosmological Perturbations. The vacuum metamorpho- 
sis scenario and its generalizations have dramatic impli- 
cations for gravitation. The transition occurs when the 
order parameter, R 2 in the TV = 4 case, reaches a crit- 



ical value. On large scales, when the cosmological cur- 
vature drops to R = m 2 the transition takes place. On 
sub-horizon scales, on the scales of voids, this critical 
value is reached a bit earlier than on average; on the 
scale of clusters, this value is reached a bit later. Since 
R = 8irG(p + 3p) in Einstein's gravitation, a high density 
or pressure will keep R above the critical value. Inside a 
cluster and on the scale of galaxies, where the mean den- 
sity is well above the cosmic background, the transition 
never takes place. 

What happens to gravitation after the transition? The 
field equations become higher-order, so that a static 
potential obeys a fourth order equation, (— AV 2 V 2 + 
V 2 )<& = AttBGp. \[A is a length, below which the poten- 
tial is exponentially suppressed. B modifies the strength 
of gravity, and grows tiny upon the transition. This 
would look like a disaster, except that the transition takes 
place only on large scales, where this Newtonian analysis 
is invalid. 

To treat the scenario in cosmological perturbation the- 
ory, we first observe that the gravitational transition 
takes place on a spacelike hypersurface of constant cur- 
vature. And the duration of the transition is effectively 
instantaneous, based on a numerical modeling of the vac- 
uum effects 6j. Taking into account small cosmological 
perturbations, this surface will be a surface of constant 
R + 5R to linear order. At times after the transition, 
the feedback of the vacuum effects on the gravitational 
field equations forces R to a constant. It is reasonable 
to surmise that fluctuations SR are forced to vanish. To 
deal with the evolution of fluctuations across the tran- 
sition we resort to junction conditions, ensuring energy 
and momentum flow is continuous. Following Ref. [3l|, 
we can choose a gauge and then match the perturbation 
variables from pre- to post-transition. In general, we can 
expect that the pure growing mode which dominates the 
evolution before the transition will give way to a linear 
combination of growing and decaying modes afterwards. 
The evolution of the perturbations after the transition, 
however, presents a challenge. 

If the transition forces fluctuations of the scalar curva- 
ture to vanish, then we find the constraint 

SR = 8nG(5p + 3Sp) 

= (h" + 3Hti - 4k 2 rj) =0 (12) 

which should be valid on the range of scales for which 
the mean curvature has frozen at R = rh 2 . Similar con- 
straints arise for the N = 2, 3 cases. Here we work in the 
synchronous gauge with metric perturbation variables h 
and 7], following the notation of Ref. H2, H = a' /a, and 
the prime indicates a derivative with respect to conformal 
time. Since density perturbations in the non-relativistic 
matter respond to gravitational fluctuations according to 
8' m = —h'/2, we have 

C + 3HS' m + 8k 2 v = 0. (13) 
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Compared to the standard case, for which 5'^ + 2Tt5' m + 
4k 2 r/ = 0, we expect stronger Hubble damping but also 
a stronger source. We are tempted to use one of the 
perturbed Einstein's equations to replace r\ with a fluid 
variable such as bp for the density perturbations. How- 
ever, we must not forget to include all contributions to 
the fluid perturbations. We know that the vacuum ef- 
fects must contribute an equivalent energy density bp v 
and pressure bp v such that Sp v + Sp m + 3bp v — 0. But 
that's all we know without making a detailed calculation. 
(A brief glance at the stress-energy tensor for the vacuum 
metamorphosis model should convince the reader that 
this is not so straightforward.) Yet, the evolution equa- 
tions for the cosmic expansion look similarly difficult and 
a simple result, transition to constant curvature, is ob- 
tained. There may be a simple resolution to the pertur- 
bation evolution as well. 

Discussion. The future prospects for distinguishing a 
gravitational transition from other proposals for dark en- 
ergy such as a cosmological constant are illustrated in 
Fig. |SJ Here, we consider the N = 4 gravitational tran- 
sition model including the contribution of scalar field 
potential energy prior to the transition. We assume a 
three-parameter family of models, consisting of Q m today 
and at the transition, and the Hubble parameter h. We 
have calculated the 95% confidence region based on the 
forecasts for CMB and SN experiments, projected into 
the fl m (transition) — fl m (today ) plane. We use the pro- 
posed SNAP experiment 33, 34] as the basis for forecasts 
to measure the recent cosmological expansion history. 
The results are shown for two different underlying mod- 
els which are currently indistinguishable, ACDM with 
ft m = 0.3 and the N = 4 model with Q m = 0.67, 0.27 
at the transition and today. In each case the Hubble pa- 
rameter has been chosen so the two models have identical 
CMB anisotropy, with the same matter density, fl m h 2 , 
and angular-diameter distance to the last scattering sur- 
face, T>cmb- Next, we expect the Planck CMB [35j ex- 
periment will determine T>cmb to 0.2% and fl m h 2 to 
0.9%, using temperature and polarization data pm 
A weaker constraint, due to our lack of knowled ge o f the 
spatial curvature, whereby a(\n VL m h 2 ) = 0.018 |38| and 
o-Qjo.'Dcmb) ~ jO-(\nQ m h 2 ) 37], is also shown. The fig- 
ure clearly indicates that these two sample cases are dis- 
tinguishable. Other probes of cosmic evolution, such as 
weak lensing and baryon acoustic oscillations can further 
sharpen the distinction. 

To summarize, we have introduced a new scenario 
which generalizes the vacuum metamorpohsis model. In 
this scenario, gravitation undergoes a transition in which 
an order parameter built out of curvature tensors freezes 
at a constant value. After the transition, the matter 
content of the universe no longer determines the cosmo- 
logical evolution — the expansion is ruled by the value 
of the order parameter. The onset of cosmic acceleration 




Q M (transition) 

FIG. 8: Projected likelihood contours are shown for Planck 
CMB and SNAP SN constraints on the N = 4 gravitational 
transition model, including the contribution of scalar field po- 
tential energy prior to the transition. Models must lie below 
and to the right of the black line, outside the gray shaded 
region, with Sl m (today) < Q m (transition). The dark red con- 
tour shows the CMB 95% likelihood region, while the blue 
line indicates the family of geometric degeneracy models. The 
lighter red contour shows a weaker constraint from the CMB 
allowing for our uncertainty in the amount of spatial cur- 
vature. The black circles show the location of two specific 
models: a ACDM model, and an N = 4 gravitational transi- 
tion. The Hubble parameters for these two models have been 
chosen so that they have identical CMB anisotropy spectra. 
The yellow contours show the SN 95% likelihood regions for 
these two underlying models. Wherease both models are vi- 
able at present, they are clearly distinguishable with future 
observations. 



is sudden, and the effective dark energy equation-of-state 
is strongly negative, w < — 1. 

We have demonstrated that a range of these mod- 
els satisfy a number of the standard tests of cosmology. 
There is excellent motivation to study these models fur- 
ther. The primary focus will be to analyze cosmological 
perturbations and the impact on structure formation. We 
can expect to find effects on the rate of growth of struc- 
ture, the large-angle CMB anisotropy pattern, and weak 
gravitational lensing. The other focus will be to exam- 
ine probes of the cosmic expansion su ggesting the dark 
energy equation-of-state sharply 0, liflkll l42l ] dropped 
below —I 0,0,0], a distinct signature of this model. 
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